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:【 目 的 ] 建 立 黑 须 污 晶 Wohlfahrtia magnifica EM ME 挖 气 黑 须 污 晶 基 因数 据 , 为 更 深 


ehe 双 峰 驼 阴 道 蝇 蛆 病 提供 理论 依据 。【 方 法 ] 采 用 高 通 量 测序 平台 lumina Hiseq™4000 对 
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黑 须 污 蝇 幼虫 \ 肾 和 成 虫 进行 转录 组 测序 ,并 进行 生物 信息 学 分 析 。 利 用 黑人 须 污 蝇 幼 虫 、 肾 和 成 忠 
录 组 数据 ,通过 基因 差异 表达 分 析 以 及 GO 功能 显著 性 富 集 分 析 的 方法 筛选 出 噢 觉 相 关 基 因 , 并 
前 过 荧光 定量 PCR 技术 对 黑人 须 污 蝇 幼虫 、 肾 和 成 蛙 中 OBP99b, OBP56a 和 OBP99a 3 个 嗅觉 相关 
因 表 达 水 平 进行 验证 。【 结果 结果 显示 ,每 个 黑 须 污 蝇 样本 的 转录 组 数据 量 在 4.96 Gb 以 上 ， 


G+C 含 量 在 35.35% 以 上 ,Q20 含量 在 97% 以 上 。 与 NCBI_nr, GO, KEGG, Pfam, Swiss-Prot 和 
eggNOG 数据 库 进行 对 比 , 共 注释 到 73 303 条 unigenes。 其 中 eggNOG 数据 库 注 释 到 的 unigenes 最 


多 ,为 35 066 条 ,分 布 在 23 个 蛋白 功能 中 ,其 中 的 翻译 修饰 、 蛋白 质 周转 的 unigenes 所 占 比 例 较 大 


GO 数据 库 注 释 到 的 unigenes 数 为 29 193 条 ,功能 包括 分 子 功能 、 细 胞 组 分 .生物 学 过 程 三 大 类 so 


个 


分 支 ,其 中 参与 生物 学 过 程 和 氧化 还 原 过 程 的 unigenes 比例 较 大 ;KEGG 数据 库 注释 到 的 


数 为 15 068 条 ,其 中 参与 信号 转 导 过 程 的 unigenes 比例 较 大 。 上 述 数 据 表 明 该 转录 组 测序 


oir 为 后 续 研 究 芮 定 了 基础 。 肾 、 成 虫 转录 组 数据 筛选 到 30 个 嗅觉 相关 
,其 中 包含 9 个 气味 结合 蛋白 (OBP) 基 因 , 进 一 步 基因 注释 发 现 OBP99b, OBP56a 和 OBP99a 


pU 污 晶 不 同 发 育 阶 段 表达 差异 显著 ( o >1, P«0.05) ,这 在 荧光 定量 PCR 分 析 结 
果 中 得 到 了 验证 。【 结 论 ] 本 研究 获得 了 黑 须 污 晶 幼虫 、. 肾 和 成 虫 转录 组 数据 ,筛选 出 黑 须 污 晶 各 
发 育 阶段 嗅觉 相关 基因 ,并 验证 了 OBP99b, OBP56a 和 OBP99a X A 4 X. 8 15 388,52 e 9 do X R3 
个 发 育 阶 段 差 异 表达 , 7] Er 26 XO JE D] :8 98 9L QE DET 5X 


X 
m 


键 词 : 黑 须 污 蝇 ; 转录 组 ; 高 通 量 测序 ; 生物 信息 学 分 析 ; 荧光 定量 PCR 
图 分 类 号 : Q966 文献 标识 码 : A 文章 编号 : 0454-6296 (2019 )06-0685-09 


Analysis of the transcriptomes and olfaction-related genes of Wohlfahrtia 


magnifica (Diptera: Sarcophagidae) at different developmental stages 
WANG Chao', BAO Hua-Er', WUYUN Gaowa’, CUI Altenwula!, BAYIN Jiragara^' , ERDEMTU" * 


(1. 


Key Laboratory of Clinical Diagnosis and Treatment Technology in Animal Disease, Ministry of 


Agriculture, College of Veterinary Medicine, Inner Mongolia Agricultural University, Hohhot 010018, 
China; 2. College of Life Science and Technology, Inner Mongolia Normal University, Hohhot 010022, 


China; 3. Barenbelle Animal Health Supervision Station of Alxa Left Banner, Bayanhot, Inner Mongolia 


75 


0308, China) 











基金 项 


: 国家 自然 科学 基金 项 目 (31360591 ) 




















作者 简介 : 





“通讯 作 
































EE, 男 1992 年 11 HÆ, 山西 长 治 人 , 硕士 研究 生 , 研究 方向 为 黑 须 污 蝇 形态 及 基因 测序 , E-mail; 779798696@ qq. com 
者 Corresponding authors, E-mail; bayin66@ 163. com; eedmt@ imau. edu. cn 





























收 稿 日 期 Received; 2018-12-21; 接受 日 期 Accepted; 2019-03-29 





686 


昆虫 学 报 Acta. Entomologica Sinica 


6235 





Abstract; [ Aim] To establish the transcriptome database of the spotted flesh fly, Wohlfahrtia magnifica , 
and to explore its genetic data, so as to provide a foundation for further studying the bactrian vaginal fly 
plague. [Methods] The transcriptomes of the larva, pupa and adult of W. magnifica were sequenced by 
a high-throughput sequencing platform (Illumina Hiseq 4000) and subjected to bioinformatics analysis. 
Olfaction-related genes were screened by differential gene expression analysis and GO function significant 
enrichment analysis based on the transcriptome data of the larva, pupa and adult of W. magnifica. The 
expression levels of three olfaction-related genes ( OBP99b, OBP56a and OBP99a) in the larva, pupa 
and adult of W. magnifica were verified by qPCR. [Results] The results showed that each sample of W. 
magnifica generated the transcriptome data of above 4. 96 Gb, in which the G + C content was above 
35. 35% , and the Q20 content was above 9796. A total of 73 303 unigenes could be annotated in NCBI. nr, 
GO, KEGG, Pfam, Swiss-Prot and eggNOG databases. Among them, the largest number of unigenes 
(35 066) were annotated in the eggNOG database, and involved in 23 protein functions, among which 
the unigenes involved in the translation modification and protein turnover accounted for a larger 
proportion. A total of 29 193 unigenes were annotated in the GO database, and involved in 50 branches 
of molecular function, cellular component and biological process, with a higher proportion of unigenes 
involved in biological process and oxidation-reduction process. There were 15 068 unigenes annotated in 
the KEGG database, with higher proportion of unigenes involved in the signal transduction process. The 
above data indicated that the transcriptome sequencing result was good, laying a foundation for the follow- 
up experiment. Based on the transcriptome data of the larva, pupa and adult of W. magnifica, 30 
olfaction-related genes including 9 odor-binding protein ( OBP) genes were screened. Further gene 
OBP99b,  OBP56a, and OBP99a 
( Ilog,FC| » 1, P «0.05) in different developmental stages of W. magnifica, which was verified by the 


annotation revealed that showed differential expression 
results of fluorescence quantitative PCR. [Conclusion] In this study, the transcriptome data of larva, 
pupa and adult of W. magnifica were obtained, the olfaction-related genes at different developmental 
stages of W. magnifica were screened, and the differential expression of OBP99b , OBP56a and OBP99a 
in the three developmental stages of larva, pupa and adult of W. magnifica was verified. These results 
provide a new idea for preventing and treating the vaginosis of the bactrian alpaca. 

Key words: Wohlfahrtia magnifica; transcriptome ; high-throughput sequencing; bioinformatics analysis; 


fluorescence quantitative PCR 


黑 须 污 蝇 Wohlfahrtia magnifica 主要 引起 母 驼 
阴道 蝇 蛆 病 , 严 重 影响 患 病 驼 的 生长 发 育 和 生产 性 
能 ,每 年 导致 阿拉 善 双 峰 驼 患 病 感染 的 发 病 率 在 
20% ~ 3096 ,死亡 率 在 2 多 左右 ( 呼 和 巴 特 尔 等 ， 
1994) 。 该 病 是 由 黑 须 污 蝇 在 骆驼 阴道 附近 排卵 ， 
在 阴道 内 寄生 并 逐渐 发 育成 1, 2 和 3 龄 幼虫 (从 1 
龄 幼虫 到 3 龄 幼虫 需 7 ~9 d) ,在 幼虫 时 期 以 母 驼 阴 











近年 来 迅速 发 展 的 高 通 量 测 序 技术 可 为 昆虫 生 
物 学 基础 研究 提供 基因 组 数据 库 , 深入 研究 发 现 生 
物 学 性 状 表达 和 行为 密切 关联 ,并 有 效 指导 得 查 发 
现 和 研究 与 目的 性 状 表达 、 遗 传 调控 相关 的 关键 基 
以 及 关键 调控 蛋白 的 功能 。 果 蝇 Drosophila z H 
前 被 人 类 研究 较为 深入 的 模式 昆虫 ,2000 年 公布 了 
全 基因 组 的 序列 框架 图 ( 刘 素 宁 和 沈 杰 , 2011), IR 






































道内 壁 为 食 对 其 造成 损伤 ,严重 的 导致 死亡 ( 额 尔 
ARESE, 2012)。 近 几 年 国外 相继 发 生 过 骆驼 感 
染 黑 须 污 晶 晶 蛆 的 病例 ,但 并 没有 找到 科学 的 治疗 
方法 (Pirali Kheirabadi et al., 2014) 。 深 入 研究 黑 须 
污 晶 的 生物 学 特征 ,通过 寻找 生物 防治 和 药物 防治 
的 新 靶 标 ,对 发 现 和 制定 该 牲畜 疾病 的 防 控 新 策略 
有 重要 的 科学 意义 。 
































觉 是 飞行 类 昆虫 克 食 .交配 \ 产 卵 行 为 的 基本 功能 之 
一 (Field et al., 2010) ,国内 外 学 者 通过 不 断 地 研究 
发 现 了 影响 果 蝇 嗅觉 发 育 的 气味 结合 蛋白 (odorant 
binding proteins, OBP) ,促进 了 对 果 晶 嗅觉 功能 的 
深入 研究 ( 吴 仲 南 等 , 2009 ) 。 根 据 果 晶 OBP 氨基 
酸 序列 的 不 同 ,将 OBP 分 为 COBP1 和 GOBP2 (Yu 
et al., 2007) ,其 中 OBP99b ( Wang et al., 2010) , 
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OBP56a ( Kumaran et al., 2018) 和 OBP99a ( Zhang et 
al., 2018) 属于 OBPs 家 族 成 员 。 对 果 晶 嗅觉 和 蛋白 
OBPs 家 族 成 员 的 深入 研究 ,也 促进 了 对 其 他 昆虫 吕 
觉 基 因 功 能 的 研究 ,例如 ,在 果 蝇 触角 中 发 现 普 通气 
味 结合 蛋白 OBPs 家 族 成 员 的 相似 序列 ,也 相继 在 
鲜 诸 日、 双 诸 日 等 昆虫 中 发 现 ( Vogt and Riddiford , 
1981 ) 。 这 些 研 究 发 现 除 OBPs 外 ,还 有 多 种 嗅觉 的 
相关 和 蛋白 参与 了 昆虫 嗅觉 识别 的 协同 功能 ,包括 气 
味 受 体 ( odorant receptors, ORs ) 和 化 学 感受 蛋白 
(chemosensory proteins, CSPs) ( 李 卫 华 等 , 2006 ) 。 
目前 国内 外 未 见 对 黑 须 污 晶 进 行 过 转录 组 测序 
分 析 及 相关 基因 功能 分 析 的 报道 。 本 研究 采用 高 通 
量 测序 技术 ,初步 建立 黑 须 污 晶 幼虫 . 肾 和 成 虫 发 育 













































































1.4 cDNA 文库 的 构建 和 Illumina 测序 

黑 须 污 晶 幼 虫 . 师 、 成 虫 cDNA 文库 的 构建 和 
Illumina 测序 由 杭州 联 川 生物 技术 股份 有 限 公 司 协 
助 完 成 。 

文库 构建 和 Inumina 测序 的 策略 是 :对 前 述 制 
备 提 取 的 黑 须 污 晶 幼虫 . 师 、 成 虫 RNA 样本 使 用 连 
接 有 Oligo ( dT) 的 磁 珠 富 集 其 mRNA ,经 抽 提 的 
mRNA 被 片段 化 试剂 随机 打 断 成 短片 段 ,以 片段 化 
的 mRNA 为 模板 ,用 六 碱 基 随 机 引物 ( random 
hexamers) 合成 cDNA 第 1 链 , 随后 加 入 缓冲 液 、 
dNTPs .RNaseH 和 DNA Polymerase I 进行 cDNA 第 
2 链 合 成 。AMPure XP beads 纯化 双 链 产物 ,利用 
T4 DNA RAHA Klenow DNA 聚合 酶 将 DNA 的 粘 


























阶段 转录 组 数据 库 , 进 而 分 析 黑 须 污 蝇 这 3 个 发 育 
阶段 与 嗅觉 功能 相关 的 3 个 基因 差异 表达 特征 ,为 
今后 深入 研究 黑 须 污 蝇 嗅觉 机 制 黄 定 基础 。 





1 材料 与 方法 


1.1 样品 的 采集 

每 年 的 5 月 底 至 10 月 初 为 阿拉 善 双 峰 驼 给 患 
阴道 蝇 蛆 病 的 多 发 时 间 段 ,随机 选 30 峰 患 病 阿拉 善 
双 峰 驼 , 在 其 阴道 处 采集 1 -3 龄 幼虫 (1 龄 第 1 -2 
天 幼虫 , 体 长 3 ~5 mm; 2 龄 第 2 -5 天 幼虫 , 体 长 
9 ~11 mm; 3 龄 第 5 -9 天 幼虫 , 体 长 13 ~ 15 mm) 
各 100 头 ,将 50 头 3 龄 幼虫 进行 实验 室 培育 ,在 3 
龄 幼虫 化 师 及 成 虫 阶 段 采 集 师 成 虫 各 20 ko ER 
采集 的 样本 需 放 入 液 氮 内 冻 存 , 放 入 液 氮 冷冻 前 需 
用 PBS 溶液 (pH 7.4, 浓 度 0.01 mol/L) 将 其 清洗 干 
净 且 不 能 有 损伤 。 
1.2. 主要 试剂 

PBS 缓冲 液 (pH 7.4, 浓 度 0.01 mol/L) (iğ 
赛 生物 公司 ) ,Dneasy Blood & Tiussue Kit (50) 试剂 
Æ ( QIAGEN Gmbh, 德国 ), 反 转录 试剂 盒 
PrimeScript"" RT Reagent Kit (DRRO37A , TaKaRa 公 
5H]), PCR 试剂 盒 Premix Taq ( DRRO92A, TaKaRa 
公司 ) ,荧光 定量 试剂 盒 SYBR 9 Premix Ex Taq™ 
(Tli RNaseH Plus) Kit (DRR420A, TaKaRa 公司 )， 
DL2000 DNA Marker (3427A , TaKaRa 公司 ) 。 
1.3 RNA 提取 

采用 TRIzol 法 对 黑 须 污 晶 幼 虫 . 晴 和 成 虫 分 别 
进行 RNA 提取 , 经 NanoDrop 检测 RNA 纯度 ， 
Qubii2. 0 定量 RNA 浓度 ,以 保证 高 质量 合格 的 样品 
进行 转录 组 测序 。 



























































性 未 端 修复 为 平 末端 ,3' 未 端 加 碱 其 A 并 加 接头 ， 
AMPureXP beads 进行 片段 选择 后 用 USER 酶 降解 
RA U 的 cDNA 第 2 链 , 进 行 PCR 扩 增 获得 最 终 测 
序 文库 。 文 库 质 检 合格 后 采用 Ilumina Hiseq""4000 
进行 测序 ,测序 读 长 为 双 端 2 x 150 bp(PE150)。 
1.5 转录 组 测序 数据 组 装 分 析 
1.5.1 Unigenes 的 注释 和 功能 分 类 : 首先 ,用 
Cutadapt 和 Perl 脚本 剔除 被 污染 适 配 性 低 、 未 确定 
的 cDNA。 然后 ,使 用 FastQC 验证 序列 质量 ,包括 
clean data 的 Q20, Q30 和 GC 含量 。 所 有 下 游 分 析 
均 基 于 高 质量 的 清洁 数据 。 用 Trinity 2.4.0 进行 转 
录 组 从 头 装配 。 将 所 有 组 装 的 Unigenes 与 非 元 余 
( Nr) 蛋白质 在 NCBI. nr, GO, KEGG, Pfam, Swiss- 
Prot 和 eggNOG 数据 库 中 进行 比 对 。 
1.5.2 测序 质 探 总 览 : 测序 产生 的 原始 数据 (raw 
data) 需要 进行 预 处 理 ,使 用 Cutadapt 去 除 测序 接头 
后 ,使 用 fqtrim 过 滤 掉 不 合格 的 序列 后 得 到 有 效 数 
据 (clean data) ,进行 下 一 步 分 析 。 

首先 去 除 测序 reads 中 的 接头 序列 ,再 对 测序 
reads 进行 窗口 法 质量 扫描 ,扫描 窗口 默认 为 6 bp, 
当 窗 口内 平均 质量 值 低 于 20 时 ,将 read 从 窗口 起 
始 到 3 ' 终 止 的 部 分 截 掉 ,去除 poly AZT ,去除 截 短 后 
长 度 小 于 100 bp 的 序列 ,去 除 截 短 后 N 的 含量 在 
5% 以 上 的 序列 。 统 计 原 始 测序 量 , 有 效 测序 量 ， 
Q20, Q30 和 GC 含量 ,并 进行 综合 评价 。 
1.6 Unigene 的 差异 表达 分 析 

用 Salmon 计算 unigenes 在 不 同时 期 的 表达 量 
(TPM) ,差异 表达 unigenes 通过 R 包 (edgeR ) 选择 
log, (fold change) > 1 或 者 log, (fold change) < -1 
且 具 有 统计 学 显著 性 (P<0.05)。 接 下 来 ,Per 脚 
本 再 次 对 差异 表达 的 unigenes 进行 CO 和 KEGG 富 
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集 分 析 ( Martin, 2011; Buchfink et al., 2015 ; Patro et 
al., 2017) 。 
1.7 荧光 定量 PCR 检测 基因 表达 

AET 1.6 节 筛 选 出 来 的 在 黑 须 污 蝇 幼虫 . 晴 、 成 
虫 中 表达 量 高 且 具 有 统计 学 显著 性 的 嗅觉 相关 基因 
OBP99b , OBP56a 和 OBP99a ,使 用 软件 Primer 5. 0 
设计 荧光 定量 引物 ( 表 1) ,引物 由 上 海 生 工 合成 ,其 
中 RPS15 为 内 参 基因 。 取 1.3 节 提 取 的 黑 须 污 晶 
ZE A RE RNA 进行 反 转录 合成 cDNA , 反 转 录 
过 程 按 照 TaKaRa 反 转 录 试 剂 盒 使 用 说 明 书 进行 。 
反应 液 在 冰 上 配制 , 反 转 录 反 应 体系 :5 x PrimeScript 
Buffer 2 (for Real Time) 4 uL, RT Primer Mix 1 uL, 















































PrimerScript RT Enzyme Mix 1 uL, Total RNA 500 
ng, 加 RNase Free dH,O 至 20 uL, PCR 程序 :37%C 
15 min, 85%C 5 s, PCR 结束 后 迅速 将 cDNA 产物 置 
F -20C 冰箱 保存 备用 。 采 用 SYBR Green 染料 法 
进行 荧光 定量 PCR, PCR 反应 体系 (20 uL): cDNA 
模板 1 pL， 上 下 游 引物 (0.4 pmol/L) 各 0.8 uL, 
SYBR 9 Greenl PCR Mix 10 uL, ddH,O 7. 4 uL, 
PCR 反应 程序 : 95C 预 变性 30 s; 95Y 5 s, 退火 温 
度 20 s, 72%C 20 s, 40 个 循环 ; 55/70 5s 289570 5s 
连续 扫描 ,0.5Y 作为 1 个 梯度 ,扫描 0.5 s。 黑 须 污 
晶 1-3 龄 幼虫 晴 、 成 虫 各 设置 5 个 生物 学 重复 ,每 
个 样本 进行 3 次 重复 实验 。 






































表 1 引物 序列 
Table 1 Primer sequences 
基因 引物 序列 (5 -3') 退火 温度 (Y ) 产物 长 度 (bp) 
Gene Primer sequences Annealing temperature Product length 
F: TTGGTTAAAGGCAGTGTCC 
Obp99b 62 152 
R: GGCGTATGGGTAAATGGTT 
F: TGCTTTCCCAACAACATGCTGC 
Obp56a 62 214 
R: TTCACCCTTGACATCCTTGC 
F: CCGCAAGGAATGTATGGAAA 
Obp99a 62 214 
R: ATCAGCCTTAACTTCGTCTT 
F: CCTACAAACCCCTTAAGCAC 
RPS15 60 224 
R: GGACCACTCGCAGAAATACA 
1.8 数据 分 析 中 注释 到 了 29 193 条 unigenes, 在 KEGG 数据 库 中 


目的 基因 的 表达 水 平 按照 内 参 基因 的 2-^"? 法 
进行 计算 , 即 ACt = Cty — Cty sag, ( Livak and 
Schmittgen, 2001) , XE [S Eoo] de 3 f JH] 2 ^" 3v, 
同一 基因 不 同时 期 相对 表达 量 方差 分 析 采 用 
SAS9.0 的 ANOVA 进行 ,Duncan 氏 法 对 平均 值 差 异 
显著 性 进行 配对 检验 。 














2 结果 


2.1 黑 须 污 蝇 幼 虫 . 肾 和 成 虫 转录 组 序列 分 析 
通过 对 黑 须 污 晶 幼 虫 . 晴 、 成 虫 转录 组 测序 , 结 
果 得 到 每 个 样本 的 碱 基数 据 量 均 在 4.96 Gb 以 上 ， 
G+C REA 35.35% 以 上 ,Q20 含量 在 97 多 以 上 。 
通过 Trinity 处 理 得 到 的 数据 ,unigenes 的 平均 长 度 
达到 了 394 nt ,表明 转录 本 建 库 ` 测 序 质 量 恨 好 ,组 
装 完整 性 较 高 ( 表 2) 。 

将 73 303 条 unigenes 5 NCBI_nr, GO, KEGG, 
Pfam, Swiss-Prot 和 eggNOG 等 6 个 数据 库 进 行 对 
比 , 共 注释 到 73 303 条 unigenes ,其 中 在 GO 数据 库 






























































得 注释 到 了 15 068 条 unigenes ,在 Pam 数据 库 中 注 
释 到 了 26 550 条 unigenes ,在 Swiss-Prot 数据 库 中 注 
释 到 了 25 040 条 unigenes ,在 eggNOG 数据 库 中 注 
释 到 了 35 066 条 unigenes, 在 NR 数据 库 中 注释 到 
T 33 399 条 unigenes( 表 3)。 

Unigene 注释 到 NR 数据 库 中 的 物种 分 布 见 图 
1 ,其 中 与 铜绿 晶 Lucilia cuprina 的 相似 基因 序列 最 
多 , 占 比 31.76% ,表明 铜绿 蝇 与 黑 须 污 晶 的 基因 具 
有 和 较 高 的 同 源 性 ,其 他 分 别 是 中 华 菊 头 蝠 
Rhinolophus sinicus (15. 66% )、 智 人 Homo sapiens 
(6.0196) , Z& W& Musca domestica ( 4. 6596 ) , JEK ok hE 
Stomoxys calcitrans (4. 6796) , “E R X ha Rhagoletis 
zephyria ( 2.71% ) 和 没有 匹配 到 同 源 性 的 物种 
(34. 54% ) 。 

与 GO 数据 库 对 比 有 3 个 本 体 (ontology ) ,分 别 
描述 基因 的 分 子 功能 (molecular function) 、 细 胞 组 分 
( cellular component ) 和 生物 学 过 程 (biological 
process) ,50 个 亚 类 ( 亚 类 是 指 可 以 在 GO 数据 库 中 
注释 到 的 蛋白 功能 ,下 同 ) ,以 这 3 种 功能 为 主体 再 





















































LW1 -3: 分 别 为 1 -3 龄 幼虫 1st -3rd instar larva, respectively; pa; "fj Pupa; ad; WIE Adult. Q20% : 
率 小 于 0.01 ) The proportion of bases with a mass value =20 (sequencing error rate is less than 0.01) ; Q30% : 








6 期 TOES. 黑 须 污 晶 不 同 发 育 阶 段 转 录 组 及 嗅觉 相关 基因 分 析 689 
m2 ， 黑 须 污 蝇 转录 组 数据 评估 
Table 2 Statistical analysis of the transcriptome data of Wohlfahrtia magnifica 
样本 原始 reads 数 ae 有 效 reads 数 We 有 效 reads 所 占 比 例 (% ) Q20 Q30 mé 
"e Number of CE Number of ME Proportion of a (33 t) 
raw reads CNN valid reads vid bae valid reads 

LW1_1 56 132 902 8.42 55 284 070 8.14 98.49 98. 72 96.50 39.20 
LW1.2 52 179 840 7.83 51 395 168 7.59 98.50 98.83 96.76 41.78 
LW1.3 42 472 184 6.37 4] 861 250 6.18 98.56 98.90 96.92 43.49 
LW1.4 44 678 044 6.70 44 104 244 6.51 98. 72 98.87 96.89 41.51 
LWI1.5 55 418 876 8.31 54 290 020 7.99 97.96 98.59 96.24 42.31 
LW2 1 46 438 314 6.97 45 799 204 6.75 98. 62 98.72 96.52 42.58 
LW2 2 42 819 500 6.42 4] 988 612 6.18 98.06 98.63 96.34 43.66 
LW2 3 46 742 926 7.01 46 146 852 6.81 98. 72 98.81 96.75 41.97 
LW2 4 44 926 360 6.74 44 062 864 6.50 98.08 98.79 96.73 42.66 
LW2_5 47 148 552 7.07 46 483 032 6.86 98. 59 98.75 96. 58 42.38 
LW3 1 45 278 048 6.79 44 706 616 6.61 98. 74 98.89 96.94 42.62 
LW3 2 58 281 970 8.74 57 639 406 8.52 98.90 98.9] 97.07 42.28 
LW3.3 47 042 638 7.06 46 383 056 6.86 98. 60 98.96 97.10 42.71 
LW3_4 38 239 182 5.74 37 815 196 5.60 98. 89 98.97 97.12 43.17 
LW3. 5 4] 783 880 6.27 4] 192 502 6.08 98.58 98. 80 96.74 43.27 
pal. 1 43 151 570 6.47 42 117 122 6.22 97.60 98.42 95.03 38.32 
pal. 2 37 763 214 5.66 36 314 672 5.37 96.16 98.50 95.13 38.33 
pal. 3 44 313 940 6.66 43 468 564 6.45 98.09 98.63 95.59 37.71 
pal_4 43 499 756 6.52 42 880 208 6.37 98.58 98.65 95.35 38.42 
pal. 5 54 789 032 8.22 46 498 688 6.76 84. 87 97.42 93.54 39.87 
ad2 1 37 671 852 5.65 37 2772 370 5.53 98.94 98.69 95.37 36.32 
ad2 2 36 100 588 5.42 35 564 930 95227 98.52 98.53 94.9 35.35 
ad2 3 34 578 160 5.19 33 461 258 4.96 96. 77 98.63 95.28 36.97 
ad2 4 42 014 678 6.30 40 231 114 5.96 95.75 98.55 95.21 37.59 
ad2_5 41 397 066 6.21 40 263 314 5.97 97.26 98. 53 95.02 36.41 


质量 值 =20 的 碱 基 所 占 比例 (测序 错误 
质量 值 >=30 的 碱 基 所 占 比 例 (测序 


错误 率 小 于 0.01) The proportion of bases with a mass value 三 30 (sequencing error rate is less than 0.01). 


A3 黑 须 污 蝇 幼 虫 . 肾 和 成 虫 转录 组 注释 结果 统计 
Table 3 Transcriptome annotation statistics of 
Wohlfahrtia magnifica larva, pupa and adult 
注释 unigenes 所 占 比例 (% ) 


注释 unigenes 数 














数据 库 l 
Number of annotated Proportion of annotated 
Database 
unigenes unigenes 
GO 29 193 39. 83 
KEGG 15 068 20.56 
Pfam 26 550 36.22 
Swiss-Prot 25 040 34.16 
eggNOG 35 066 47.84 
NR 33 399 45.56 
总 计 Total 73 303 100 


次 对 unigene 进行 分 类 得 到 以 下 3 组 数据 (图 2) : 生 
物 学 过 程 中 含有 43 232 条 unigene,25 个 亚 类 ,其 中 


生物 学 过 程 和 氧化 还 原 过 程 所 占 的 比例 最 大 ,unigene 
数 达 到 了 1 500 以 上 ;细胞 组 分 中 含有 21 920 条 
unigenes,15 个 亚 类 ,其 中 细胞 核 和 细胞 质 的 unigenes 
占 比 例 最 大 ,unigene 数 超过 4 000; 分 子 功能 中 含有 
21 721 条 unigenes ,10 个 亚 类 ,其 中 蛋白 结合 unigenes 
所 占 比例 最 大 ,unigene 数 超过 了 3 000。 根 据 以 上 结 
果 得 出 在 黑 须 污 昌 幼虫 、 肾 成 虫 发 育 过 程 细 胞 中 化 
学 反应 过 程 以 及 和 蛋白质 转录 翻译 的 重要 性 。 

KEGG 富 集 分 析 表 明 , 黑 须 污 蝇 幼虫 、 肾 和 成 虫 
unigenes 富 集 在 了 6 个 大 类 和 40 个 亚 类 ( 亚 类 是 指 
可 以 在 KEGG 数据 库 中 注释 到 的 具体 的 信号 通路 ) 
(图 3) ,其 中 信号 转 导 富 集 的 unigenes 条 数 最 多 , 达 
到 了 2 164 条 。 说 明 信 和 号 转 导 在 黑 须 污 蝇 幼虫 . 师 、 
成 虫 发 育 中 有 重要 的 作用 。 
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1 黑 须 污 晶 幼 虫 师 、 成 虫 unigenes 在 NR 数据 库 中 的 物种 分 布 


Species distribution of unigenes of Wohlfahrtia magnifica larva, pupa and adult in the NR database 
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2 黑 须 污 晶 幼虫 . 师 和 成 虫 unigenes GO 功能 注释 
Fig. 2 GO function annotation of unigenes of Wohlfahrtia magnifica larva, pupa and adult 





与 eggNOG 数据 库 进 行 对 比 和 假定 蛋白 功能 分 
析 ,共有 35 066 条 unigenes 被 注释 到 eggNOG 的 23 
个 蛋白 功能 条 目下 (蛋白 功能 条 目 包括 未 知 蛋 白 功 
能 可 以 和 其 他 已 知 蛋白 功能 条 目 并 列 ) (图 4)。 其 
中 共有 15 755 条 unigene 为 功能 未 知 ,说 明 有 未 知 
基因 不 能 在 eggNOG 数据 库 中 得 到 注释 ,其 次 依次 
是 翻译 后 修饰 (3 362 条 unigenes) 分泌 和 圳 泡 运 输 
(1 864 条 unigenes) .翻译 和 核糖 体 结构 (1 782 条 
unigenes) 。 表 明了 翻译 修饰 功能 、 蛋 白质 周转 在 黑 


须 污 晶 幼 虫 肾 、 成 虫 生长 发 育 中 有 一 定 的 作用 。 
2.2 在 黑 须 污 蝇 幼 虫 . 晴 和 成 虫 中 嗅觉 相关 基因 的 
差异 表达 与 GO 富 集 功能 分 析 
es E 
基因 差异 表达 分 析 以 及 GO 功能 显著 性 富 集 分 析 发 
Rod ui 
CSP41B, cspA, CspgA), OR 基因 10 个 ( OR3, 
ORIFI2, Or30a, Or63a, Orcl, Orc2, Orcà, Orc4, 
Orc5, Orc6) , IR 基因 5 ^F (125a, ird5, IRAKI, 
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黑 须 污 晶 不 同 发 育 阶段 转录 组 及 嗅觉 相关 基因 分 析 691 
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Fig. 3 Classification of KEGG pathways of unigenes of Wohlfahrtia magnifica larva, pupa and adult 
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类 别 Category 


m A: RNA 加 工 和 修饰 RNA processing and modification 

=B: 染色 质 结 和 动力 学 Chromatin structure and dynamics 
mC: 能 源 生 产 和 转换 Energy production and conversion 

"D: 细胞 周期 控制 和 细胞 分 裂 Cell cycle control and cell division 
m E: 氨基 酸 转运 和 代谢 Amino acid transport and metabolism. 

m F: 核 昔 酸 转运 和 代谢 Nucleotide transport and metabolism 

m G: 碳水 化 合 物 的 运输 代谢 Carbohydrate transport and metabolism 
mH: 辅酶 转运 和 代谢 Coenzyme transport and metabolism 

5 p: 脂 质 运 输 和 新 陈 代谢 Lipid transport and metabolism 

my 翻译 和 核糖 体 结构 Translation and ribosomal structure 

m K: 转录 Transcription 

mL: 复制 、 重 组 和 修复 Replication, recombination and repair 
mM: 细胞 壁 / 膜 / 包 膜 生 物 学 过 程 Cell wall/membrane/envelope biogenesis 
=N: 细胞 运动 Cell motility 

m O: 翻译 后 修饰 Posttranslational modification 

m P: 无 机 离子 转运 和 代谢 Inorganic ion transport and metabolism 
9 Q: 次 生 代谢 产物 的 生物 合成 Secondary metabolites biosynthesis 
m S: 功能 未 知 Function unkown 

mT 信号 转 导 机 制 Signal transduction mechanisms 

mU: 分 泌 和 襄 泡 运输 Secretion and vesicular transport 

= v: 防御 机 制 Defense mechanisms 

m y: 核 结构 Nuclear structure 

m Z: 骨架 Cytoskeleton 














4 黑 须 污 晶 幼虫 . 师 和 成 虫 unigenes 的 eggNOG 功能 分 类 





Fig. 4 


IREH, Irbp) ,SNMP 基因 2 个 (Snmpl,， Snmp2) , OBP 
基 9 个 (Obpl9a, Obp19d, OBP56a, OBP99b, 
OBP99a , Obp56d, Obp69a, Obp83b, Obp70) , diE— 2E 
基因 差异 表达 分 析 表 明 , OBP99b 基因 的 1log,FC| > 
1 且 P<0.05 ,说明 OBP99b TEE E Ts We A [n] A 3 DT 
段 的 表达 存在 显著 差异 ,通过 GO 富 集 分 析 发 现 
OBP99b 富 集 在 嗅觉 行为 和 对 信息 素 的 反应 2 个 
GO term F; OBP56a 基因 的 llogyFC1>1 且 P< 
0. 05 ,说 明 OBP56a 在 黑 须 污 蝇 不 同 发 育 阶段 的 表 
达 存 在 显著 差异 ,通过 GO 富 集 分 析 发 现 OBP56a 
富 集 在 化 学 刺激 的 感官 知觉 和 和 气味 结合 、 嗅 觉 行 为 





eggNOG functional classification of unigenes of Wohlfahrtia magnifica larva, pupa and adult 


3 个 GO-term F; OBP99a 基因 的 llogyFC1 »1 H 
P<0.05 ,说 明 OBP99a 在 黑 须 污 晶 不 同 发 育 阶段 的 
表达 存在 显著 差异 ,通过 GO 富 集 分 析 发 现 0BP99a 
富 集 在 嗅觉 行为 、 化 学 刺激 的 感官 知觉 .气味 结合 
对 信息 素 的 反应 4 个 COterm 下 。 
2.3 ”荧光 定量 PCR 结果 检测 3 个 嗅觉 相关 基因 
mRNA 表达 量 
以 不 同时 期 黑 须 污 蝇 样品 cDNA 为 模板 ,分 别 

对 待 测 基因 进行 荧光 定量 检测 。0BP99b,，0BP56a 
和 OBP99a 融 解 曲线 结果 表明 ,峰值 单一 ,各 样品 熔 
解 点 温度 相同 ,说 明 扩 增产 物 特异 性 较 强 ,反应 条 件 








692 昆虫 学 报 Acta Entomologica. Sinica 62 35 





好 ,无 引物 二 聚 体 及 其 他 非特 异性 扩 增 产物 。 奖 光 
定量 PCR 结果 表明 ，OBP992 在 黑 须 污 蝇 晴 中 表达 
量 极 显 著 高 于 幼虫 和 成 虫 中 的 (P <0.01) , OBP56a 
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在 2 龄 幼虫 中 的 表达 量 极 显 著 高 于 1 和 3 龄 幼虫 以 
及 晴 和 成 虫 中 的 (P<0.01) ,OBP99a 在 师 和 成 虫 中 
的 表达 量 极 显著 高 于 幼虫 中 的 (P<0.01) (图 5)。 
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图 5 OBP99b (A), OBP56a ( B) fll OBP99a 〈(C) 基 因 在 黑 须 污 蝇 不 同 发 育 阶 段 的 表达 量 
Fig. 5 Relative expression levels of OBP99b ( A), OBP56a (B) and OBP99a (C) at different 
developmental stages of Wohlfahrtia magnifica 
LWl, LW2, LW3: 分 别 为 1-3 龄 幼虫 1st — 3rd instar larva, respectively; pal ; Hý Pupa; ad2: 成 虫 Adult 组 内 生物 学 重复 n=5, 柱 上 不 同 大 写 
字母 表示 基因 表达 量 在 不 同 发 育 阶 段 间 差异 极 显 著 (P<0.01, Duncan 氏 检 验 ) 。Intra-group biological repeats n — 5. Different capital letters above 


bars show extremely significant difference in gene expression level among different developmental stages ( P «0.01, Duncan' s test). 








3 讨论 


本 研究 首次 对 黑 须 污 晶 幼 虫 . 晴 和 成 虫 进 行 了 
转录 组 测序 以 及 序列 组 装 与 分 析 , 共 获得 了 162. 04 
Gb 的 有 效 碱 基 , 73 303 条 unigenes。 所 有 样本 的 
Q20 和 Q30 碱 基 比 例 不 低 于 93% ,与 NR 数据 库 对 
比 得 知 铜绿 蝇 ( Anstead et al., 2015) 与 黑 须 污 晶 基 
因 具 有 较 高 的 同 源 性 ,其 中 仍 有 34.54% 的 unigenes 
没有 得 到 注释 ,可 能 是 未 知 的 基因 或 非 编 码 序列 
( Hou et al., 2011), 

GO 富 集 分 析 对 unigenes 进行 功能 分 类 有 利于 
全 面 了 解 黑 须 污 晶 幼虫 . 师 和 成 虫 的 基因 和 基因 产 
物 。 在 GO 数据 库 中 一 共 注 释 到 29 193 条 基因 ,并 
将 其 分 为 分 子 功 能 、 细 胞 组 分 和 生物 学 过 程 三 大 组 ， 
其 中 参与 生物 学 过 程 、 细 胞 核 和 细胞 质 、 蛋 白质 结合 



































的 基因 功能 注释 较 多 。 黑 须 污 晶 幼 虫 . 晴 和 成 虫 转 
录 组 经 KEGG 数据 库 注释 得 到 200 条 代谢 通路 , 主 
要 有 信和 号 转 导 等 ,说 明 信 和 号 转 导 在 黑 须 污 晶 嗅 觉 行 
为 活动 中 的 重要 性 。 

本 研究 通过 分 析 黑 须 污 蝇 幼虫 、 肾 和 成 虫 转录 
组 数据 得 知 参 与 环境 信息 处 理 中 信号 转 导 基 因数 量 
最 多 ,有 2 164 条 ,发 掘 出 30 个 嗅觉 相关 的 基因 , 包 
$54 ^ CSP 基因 ,10 个 OR 基因 ,5 个 IR 基因 2 个 
SNMP 基因 和 9 个 OBP 基因 ,充分 说 明了 黑 须 污 蝇 
各 发 育 阶 段 嗅觉 与 周围 环境 交互 的 重要 作用 。 
OBPs 与 气味 分 子 结合 是 昆虫 于 外 界 进 行 信号 传递 
的 一 个 关键 步骤 并 在 嗅觉 发 育 中 发 挥 着 重要 作用 
( Christiansen et al., 2012) ,本 研究 通过 基因 差异 表 
达 分 析 以 及 GO 富 集 分 析 发 现 9 个 OBP 基因 ,其 中 
OBP99b, OBP56a 和 OBP99a 在 黑 须 污 晶 的 3 个 发 
育 阶段 有 显著 差异 表达 (P <0.05), 并 且 表 达 量 较 
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大 ,GO 富 集 分 析 得 知 它们 全 部 富 集 在 嗅觉 行为 下 。 
最 后 经 过 荧光 定量 PCR 实验 验证 了 上 述 3 个 基因 
在 黑 须 污 晶 幼虫 . 师 和 成 虫 阶段 表达 水 平 差 异 极 显 
著 (P<0.01) ,由 此 说 明 上 述 3 个 基因 与 黑 须 污 晶 
嗅觉 的 发 育 有 关联 。 与 其 他 的 昆虫 相 比 (Crava et 
al., 2019) , 黑 须 污 晶 幼虫 中 存在 差异 表达 的 嗅觉 相 
关 基 因数 量 相对 较 少 ,这 可 能 与 其 幼虫 阶段 没有 触 
角 器 官 相 关 。 有 关 OBP 基因 在 黑 须 污 蝇 幼 虫 . 晴 和 
成 虫 中 的 特异 性 及 其 他 功能 有 待 于 进一步 的 研究 。 

本 研究 填补 了 黑 须 污 晶 转录 组 数据 库 空 白 , 验 
证 了 OBP99b, OBP56a 和 OBP99a 基因 与 黑 须 污 晶 
幼虫 . 师 和 成 虫 嗅觉 发 育 相 关联 ,并 为 基因 工程 技术 
防治 黑 须 污 晶 危 害 提 供 新 理念 。 
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